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We consider two species of bosons in one dimension near the Tonks-Girardeau limit of infinite 
interactions. For the case of equal masses and equal intraspecies interactions, the system can be 
mapped to a S' = 1/2 XXZ Heisenberg spin chain, thus allowing one to access different magnetic 
phases. Using a powerful ansatz developed for the two-component Fermi system, we elucidate the 
evolution from few to many particles for the experimentally relevant case of an external harmonic 
confinement. In the few-body limit, we already find clear evidence of both ferromagnetic and 
antiferromagnetic spin correlations as the ratio of intraspecies and interspecies interactions is varied. 
Furthermore, we observe the rapid emergence of symmetry-broken magnetic ground states as the 
particle number is increased. We therefore demonstrate that systems containing only a few bosons 
are an ideal setting in which to realize the highly sought-after itinerant ferromagnetic phase. 


Quantum magnetism is ubiquitous in nature and plays 
a central role in important phenomena such as high- 
temperature superconductivity [1]. Furthermore, it un¬ 
derpins the technological advances in data storage [2], 
and it promises a new generation of spintronic devices, 
where the spin of the electron rather than its charge is 
used to transfer information. However, despite its ubiq¬ 
uity, magnetic phenomena are often difficult to charac¬ 
terise and treat theoretically; for instance, itinerant fer¬ 
romagnetism of delocalized fermions requires strong in¬ 
teractions and is thus not completely understood [3]. 

One can gain insight into magnetic phases by con¬ 
sidering cleaner, more idealised versions of the phenom¬ 
ena. In particular, an important question is whether or 
not ferromagnetism can exist without an underlying lat¬ 
tice [4] . This possibility was investigated experimentally 
in atomic Fermi gases [5, 6], but the system proved to 
be unstable towards fermion pairing rather than mag¬ 
netism. More recently, it has been proposed that itiner¬ 
ant ferromagnetism can be realised in a strongly interact¬ 
ing one-dimensional (ID) Fermi gas [7]. However, in this 
case, one cannot access the ferromagnetic phase without 
explicitly breaking the spin symmetry with an external 
field [8] , and this potentially complicates the observation 
of ferromagnetism. 

In this Letter, we show that both itinerant ferromag¬ 
netism and antiferromagnetism can be investigated with 
a ID two-component mixture of bosons [9]. For the 
case of equal masses and equal intraspecies interactions 
(fftt = Oil), the Bose-Bose mixture may be regarded as a 
pseudo-spin S' = 1/2 system, and it can be mapped to a 
ID spin chain in the limit of strong interactions [10, 11]. 
In particular, for infinite intraspecies interactions, the 
system is formally equivalent to a ID Fermi gas. How¬ 
ever, in contrast to the Fermi case, one can break the 
SU(2) symmetry and access different magnetic phases 
by varying the ratio of intraspecies and interspecies in¬ 
teractions (see Fig. 1). Moreover, the ferromagnetic state 




FIG. 1. (a) Magnetic phase diagram for the two-component 
Bose gas near the Tonks-Girardeau limit of infinite interac¬ 
tions, llgvv' 0, where we focus on the “upper-branch” 
ground state and ignore any two-body bound states, (b) The 
system may be mapped onto an XXZ Heisenberg spin model, 
Eq. (3), with couplings and in the x-y plane and 2 di¬ 
rection, respectively. Equivalent points in the phase diagrams 
are marked by identical symbols, e.g., /> = 37r/4 corresponds 
to identical inter- and intraspecies interactions. The thick 
(red) lines correspond to the XXX Heisenberg model that is 
realized by the two-component Eermi gas. In the thermody¬ 
namic limit, these lines bound regions of Ising ferromagnetism 
(FM) and Ising antiferromagnetism (AFM), separated by dis¬ 
ordered “spin liquid” phases. In the hatched region all inter¬ 
actions are repulsive, and the spin-chain ground state is the 
actual ground state of the system. 

should be stable, unlike in higher dimensions [12-14]. 

Experimentally, ID quantum gases have been success¬ 
fully realized with strongly interacting bosons [15, 16], 
two species of fermions [17-20], and, more recently, 
fermions with SU(n) symmetry, where the number of 
spin components can be tuned from n = 2 to 6 [21]. All 
these ID experiments feature an underlying harmonic 
trapping potential, which has enabled the study of the 
evolution from few to many particles [18, 20, 22]. How¬ 
ever, it is theoretically challenging to treat ID interacting 
particles in a harmonic trap, since the problem cannot 
be solved exactly in general [23]. Only a small num¬ 
ber of trapped particles can be treated exactly numeri- 
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cally [22, 24-32], 

Recently, we developed a highly accurate ansatz for the 
strongly interacting ID Fermi gas in a harmonic poten¬ 
tial that allowed us to analytically determine the wave- 
function for a single spin-J, impurity with any number of 
majority t fermions [33], Here, we demonstrate that this 
approach can be extended to provide an extremely accu¬ 
rate description of the Bose mixture, providing a bench¬ 
mark for other numerical techniques. We use this method 
to address the question of how the ground state of the 
few-body system is connected to magnetic phases in the 
many-body limit. As the particle number N is increased, 
we observe the rapid emergence of a doubly degenerate 
ground state, corresponding to the Z 2 symmetry that is 
ultimately broken within the Ising magnetic phases in the 
thermodynamic limit. We furthermore argue that stable 
ferromagnetism is best observed in a ID few-body system. 

Effective spin model .— We start by considering a ID 
two-component {v =t,i) gas of bosons with mass m, 
trapped in a harmonic potential with frequency uj. Inter- 
and intraspecies interactions are assumed to be short- 
ranged, and parametrized respectively by g^i and g^,^. 
The Hamiltonian is 


9^ 1 2 2 I / \ 


+ I dx (1) 


where V'J and respectively, create and annihilate a bo¬ 
son with pseudo spin v. In the following we use harmonic 
oscillator units, where h = m = lo = 1. Moreover, we 
consider the symmetric case where the two intraspecies 
interactions are of equal strength and the number of par¬ 
ticles in each spin state is the same, i.e., g^^ = and 
N^ = N^ = N/2. 

In the Tonks-Girardeau limit \gu,y' \ —> 00 , the particles 
become impenetrable and retain their particular order¬ 
ing. Thus, even in the absence of an underlying lattice, 
one can consider the system as a discrete chain of length 
N, where the site index defines the position along the 
ordered chain. Perturbing away from this limit allows 
nearest neighbors to exchange position and one can de¬ 
scribe the system using a spin-chain model [11, 28, 29, 33- 
35]. To make this connection evident, it is convenient to 
replace the inter- and intraspecies interactions g^^ and 
3TT = 9ll effective parameters g > 0 and (f, such 
that the entire phase diagram is traversed by chang¬ 
ing the “angle” (/>, as depicted in Fig. 1. These are re¬ 
lated to the original interactions by siniji = g/g^\^ and 
coscf) = g/g^i — 2g/g^^. To linear order in 1/g, the origi¬ 
nal Hamiltonian may then be written as [35, 36]: 

H = ee> + [/^e + n)/g. ( 2 ) 


Here, eg is the energy of JV identical fermions, and ?{ is 


an effective Heisenberg XXZ Hamiltonian: 

N-l 

+ crlcrf+i) + JVfcrf+i] , (3) 

where J-*- = sintj), = cost/), and cr^ are the Pauli 
spin matrices at position i. The constant Ae = (J^ — 
Vi ensures that only identical (distinguish¬ 
able) bosons interact via g^^ (fftf)- ^^ote that, apart from 
parity, the properties of the eigenstates of H are invariant 
under the transformation (j) 1 —> —(f [37]. 

The energy cost of exchanging opposite spins at po¬ 
sitions i and i -|- 1 depends on the underlying external 
confinement, an effect encapsulated in the coefficients rji 
of Eq. (3). These must in general be calculated numer¬ 
ically, but a key simplification is that they are indepen¬ 
dent of the interactions and thus one need only determine 
them for a single (j). In particular, for infinite intraspecies 
interaction {(j> = tt/A, 57r/4), the system is equivalent 
to a two-component Fermi gas, where the rji have al¬ 
ready been investigated for A < 9 [28, 29, 33]. We have 
recently shown for the harmonically trapped Fermi gas 
that gi is well approximated by the analytic expression 
i{N — i)fj{N) [33]. Referring to the Supplemental Mate¬ 
rial [37], we find that this approach also yields near exact 
wavefunctions for any (j). We will, thus, use this ansatz 
to consider larger numbers of particles. 

Phase diagram. — To gain insight into the behavior of 
the trapped boson system, it is instructive to first con¬ 
sider the uniform XXZ model {gi = 1), where the mag¬ 
netic phases are well characterized [38, 39]. In the ther¬ 
modynamic limit, the local density approximation (LDA) 
shows that the topology of the uniform phase diagram is 
the same as that of the trapped system. As displayed 
in Fig. 1, we expect to have fully saturated ferromag¬ 
netism (FM) of the Ising type when i.e., 

when 37r/4 < (/) < 57r/4 [38]. Indeed, the first-order tran¬ 
sition at = — J-*- < 0 corresponds to g^^ = g^^ > 0, 
the criterion for phase separation in a weakly interacting 
two-component Bose gas [39]. However, as we discuss 
below, the interface between f and f bosons is sharper 
in the strong-coupling limit, being of order the inter¬ 
particle spacing. The other first-order FM transition, 
= J-*- < 0, has no weak-coupling analogue since it 
occurs on the metastable “upper branch” for attractive 
interactions (or “super-Tonks” regime). 

When I J-*-| = I J*|, % becomes equivalent to the XXX 
Heisenberg model that can be realized with fermions, e.g., 
(/) = 37r/4 can be mapped to the upper branch of the at¬ 
tractive Fermi gas [28, 35], which is the regime where 
itinerant ferromagnetism exists [7]. Once ^fg^v < 0, 
we will also have Ising antiferromagnetism (AFM) for 
—-KjA < (j) < 7r/4, but the phase transition to this state 
is of the Kosterlitz-Thouless type in the uniform case [38] . 

Ground-state energy. — We now turn to the question 
of how the few-body system approaches the many-body 
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FIG. 2. Shift in the ground-state energy (see main text) as a 
function of ij> for increasing particle number, from = 2 to 16 
(solid lines). The dash-dotted black line is an extrapolation 
to the limit N ^ oo [40]. The black symbols and shaded 
regions correspond to the interaction strengths and magnetic 
phases illustrated in Fig. 1. Diamonds indicate analytical 
Bethe ansatz results within LDA, where <() = 7r/4 corresponds 
to a two-component Fermi gas (4), and (j> = 37r/4 is equivalent 
to identical bosons (5). The dotted green line corresponds to 
the result for spin-polarized bosons for A —^ oo. 

limit. We first consider the ground-state energy of the 
trapped gas: Eq = Eq — C'gs/ff [41]. Since £o is a simple 
energy shift, we take Cgs as a function of (p, as shown in 
Fig. 2. The result for = 1 exhibits disconti¬ 

nuities in the slope at <() = 0, tt, corresponding to level 
crossings in the energy spectrum, where the total spin of 
the ground state is S' = 0 for 0 < (^ < tt and S = 1 for 
—TT < (j) < 0. In this case, there is no intraspecies inter¬ 
action to couple states of different total spin, similarly to 
the two-component Fermi gas, and one cannot adiabati- 
cally evolve the ground state. However, for N > 2, the 
total spin is no longer conserved and thus one can access 
different magnetic phases by varying p adiabatically. 

Figure 2 shows that the results rapidly approach a uni¬ 
versal curve with increasing N. In particular, we see 
that certain (p extrapolate to known Bethe ansatz results 
within LDA [42]. For p = 7r/4, where our system cor¬ 
responds to a two-component Fermi gas, we recover the 
LDA ground-state energy of the repulsive Gaudin-Yang 
model [42-44]: 

Egy I 128^2 ln(2) ViV 
N'^ 2 457r2 

For p = 37r/4, we have g^i, = and the ground state is 
equivalent to that of a strongly repulsive gas of N iden¬ 
tical bosons. Thus, we see that our results converge to 
the LDA ground-state energy of the Lieb-Liniger model 
at strong coupling [42]: 

£;ll 1 128^2 y/N 

2 4571^ g^'f- 


Within the entire region 37r/4 < p < 57r/4, we find 
that the ground-state energy in the limit iV —)■ oo corre¬ 
sponds to that of spin-polarized bosons, where it follows 
from Eq. (5) that Cgs = ^^ 4 ^ 2 '' cos{p — 37r/4). For fixed 
spin populations = N^, this result physically corre¬ 
sponds to phase separation of f and P particles, where we 
can neglect the energy of the interface in the thermody¬ 
namic limit since it scales like 1/N compared to the total 
energy Eq. The formation of spin-polarized domains is a 
signature of ferromagnetism [45], and the emergence of 
kinks in the energy at ^ = 37r/4 and 57r/4 clearly indi¬ 
cates a first-order transition to the FM phase (Fig. 2). 
While there are no sharp features in the energy at the 
AFM phase boundaries, the fact that Cgs = 0 at (^ = 0 
is a direct consequence of classical Neel ordering, where 
identical spins are staggered and thus experience no in¬ 
traspecies interaction. 

Magnetic phases. — A clear signature of magnetic order 
can be extracted from the degeneracy of the ground state. 
A two-fold degenerate ground state is a necessary condi¬ 
tion for Ising magnetic order, and thus its emergence in 
the few-body system is a precursor for symmetry-broken 
magnetic states as N ^ 00 . As shown in Fig. 3(a), the 
two lowest energy states become rapidly degenerate with 
increasing N in the FM regime. In the few-body sys¬ 
tem, these two states have opposite parity, but as they 
become degenerate they can be combined to form states 
that break spin symmetry and parity. Physically, this 
degeneracy emerges once the tunnelling between differ¬ 
ent symmetry broken states is suppressed for large N. 
A similar situation occurs in the AFM region, but the 
degeneracy manifests more slowly. Therefore, very few 
particles are necessary to observe strong magnetic corre¬ 
lations and a large breaking of the Z 2 symmetry in the 
FM region, while larger systems are needed to observe 
similar effects in the AFM region. 

We can further characterize the different phases using 
the short-range correlation function Xi = in the 

ground state. Referring to the insets of Fig. 3(b), we 
see that Xi is always close to —1 around p = 0, signifying 
AFM correlations. Likewise, for the FM phase, we obtain 
Xi « 1 throughout the spin chain, excluding the interface 
between the spin-polarized domains. For strong interac¬ 
tions, this domain wall is sharp, decaying exponentially 
with distance even in the few-body system. 

To obtain a more global description of magnetic cor¬ 
relations, we take the trap-averaged function x{4>) = 
, normalized such that x(0) = — 1 and x(7r) = 1. 
At the FM boundaries, the ground state is a Heisen¬ 
berg ferromagnet with total spin S = N/2 oriented in 
the x-y plane. Here, the wavefunction is obtained by 
applying the total spin lowering operator N/2 times to 
the spin-polarized system jff ... f). Thus, by symmetry, 
X = Xi = 0- With increasing N, we find that the av¬ 
eraged short-range correlations approach a rectangular 
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FIG. 3. (a) Energy-difference of ground and first excited 

states, El — Eo, compared to the energy of the second ex¬ 
cited state, i? 2 . (b) Trap-averaged magnetic correlation func¬ 
tion evaluated in the ground state. The dash-dotted line is 
an extrapolation to the thermodynamic limit [40], and the 
stars are analytic results for a uniform system. The insets 
show the spatially-resolved correlation function for = 16 
and two different values of 4>. (c) Spin gap AS for the spin 
polarised (dashed) and unpolarised (solid) ground states. The 
green dotted lines are the analytical results for large N — see 
Supplemental Material [37]. 


function in the FM region, where 1 — x — (tt — 4'Y/N 
around (^ = tt for ^ 1. This is consistent with the 
appearance of classical Ising FM along the z direction. 

By contrast, we see that quantum fluctuations reduce 
the Ising AFM correlations around (j) = 0, with 1+x — 
in the limit N —>■ oo. Furthermore, x evolves smoothly 
with (p from the AFM phase to the disordered “spin liq¬ 
uid” occupying the region 7r/4 < \(j)\ < Sir/i. The dis¬ 
ordered liquid also features x < 0 since its spin corre¬ 


lations resemble those in a Fermi system. Indeed, for 
(j) = -Ett 12, the spin model can be formally mapped onto 
free fermions using the Jordan-Wigner transformation 
[39]. For the uniform case, the free-fermion problem can 
be solved analytically and we find x = ~4/7r^, which 
compares well with the large N limit of the trapped sys¬ 
tem [see Fig. 3(b)]. 

Spin imbalance .— When ^ N\^, magnetic corre¬ 
lations strongly influence the behavior of the trapped 
system. In the limit of a single impurity = 1, we 
find that the type of correlations determines the posi¬ 
tion of the impurity: in the AFM regime, the impu¬ 
rity sits at the trap center [33], while FM correlations 
confine it to the edge. More generally, we find in the 
few-body system that the fully spin-polarised state al¬ 
ways has the lowest energy within the FM regime, as 
expected, while the ground state for all other p is un¬ 
polarised. In practice, one can observe this by adding a 
small spin-interconversion coupling [36] . To quantify this 
behavior, we determine the energy cost of flipping a spin 
in the ground state [Fig. 3(c)]. For the fully polarized FM 
ground state, we see that the “spin gap” AE always van¬ 
ishes at the FM phase boundaries and rapidly converges 
to a universal curve with increasing N [37]. This further 
supports the idea that FM can be realized in a few-boson 
system. On the other hand, for the unpolarized system 
we find that the AFM spin gap remains finite even in the 
disordered phase, illustrating the different nature of the 
spin excitations. In both cases, we find that the behavior 
in the trap is qualitatively different from that in the uni¬ 
form system, since the flipped spins predominantly reside 
at the trap edges and, hence, the spin gaps may not be 
obtained via the LDA. 

Concluding remarks .— The physics described in this 
Letter may be realized in a two-component Bose gas close 
to overlapping ID inter- and intraspecies resonances, 
which obviously necessitates some fine tuning. However, 
we emphasize that the ID Bose gas is substantially more 
tunable than its 3D counterpart: Firstly, the magnetic 
field need not be set exactly to any of the 3D resonances, 
since part of the fine tuning can be done by varying the 
strength of the transverse confinement [46] , and this may 
even be done in a species-selective manner. Secondly, 
we propose to employ any two hyperfine spin states, as 
spin-changing collisions will be suppressed in the strongly 
interacting “fermionized” limit. Thus we anticipate that 
the strongly-coupled two-component Bose gas could be 
investigated in current experimental setups. 

We have argued that clear signatures of a first-order 
quantum phase transition to a ferromagnetic state - the 
analogue of itinerant electron Stoner ferromagnetism - 
appear already in the few-body limit. Furthermore, we 
emphasize that the few-body scenario is likely to be ideal 
for observing magnetic phases since the condition of stay¬ 
ing in the strongly interacting regime, g/y/N 3> 1, be¬ 
comes increasingly difficult to satisfy as N increases. 
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Smaller values of N also enlarge the relative energy 
spacing between eigenstates, thus making it easier to 
tune interactions adiabatically and access stable mag¬ 
netic ground states that are robust with respect to ther¬ 
mal fluctuations. Hence, the harmonically trapped and 
strongly interacting few-body Bose gas appears ideally 
suited for the realisation of quantum magnetic phases. 
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SYMMETRIES OE THE HAMILTONIAN 

As discussed in the manuscript, when we are in the Tonks-Girardeau limit | —)■ cx), the system may be described 

in terms of an effective XXZ Hamiltonian, 


N-l 


n = 


i=l 


Vi («+i + CTfcjf+i) + J"cr^iCr^+i] , 


(SI) 


where the parameters are defined as in the main text. For even particle number N, as is the case when the spin 
populations are balanced {N^ = N^), one can show that this spin-chain Hamiltonian satisfies 'H(—</>) = U'H{(j))U\ 
with the unitary operator U = ^ii- Thus, with the exception of parity, we expect the properties of the eigenstates 

to be invariant under (j) !->• —cj). The parity operator P = so that for even N/2 we have [U,P] = 0 and 

the ground state always has the same parity. For odd N/2, instead {U,P} = 0, and consequently the ground state 
changes parity at 0 = 0 and tt. On the other hand, the total energy is not invariant under the unitary transformation 
due to the presence of the constant Ae{(j)) in the strong-coupling Hamiltonian, Eq. (2). 


PARAMETERS OE THE MODEL 


The parameters rji which yield the values of the spin-exchange coefficients of the spin-chain Hamiltonian (SI) were 
computed in Refs. [28, 29, 33]. For completeness, we list these in Table SI up to A = 9. 


A 

2 

3 

4 

5 

6 

7 

8 

9 

Vi 

^1/277 

27/(16^^) 

0.893823 

1.08303 

1.25109 

1.40370 

1.54440 

1.67556 

V2 



1.17325 

1.58860 

1.95105 

2.27645 

2.57407 

2.84986 

773 





2.17856 

2.70016 

3.17251 

3.60715 

V4, 







3.36929 

3.97952 


TABLE SI. Distinct parameters of the spin-chain Hamiltonian. The full set can be obtained using Tjjv-i = ?7i. For N < 3 
the coefficients are exact. For N = 4 the coefficients may be obtained analytically [33], however the resulting expressions are 
cumbersome and are not repeated here. 


In the manuscript, instead of using the numerically exact coefficients in Table SI we use an approximate model 
derived in Ref. [33], where r]i = i{N — i)fj{N). Note that fj{N) only sets the absolute value of the energy and does not 
affect the eigenstates or the energy ratios for each N. We fix 77 (A) as following: At cj) = tt, corresponding to infinite 
interspecies interactions, we know the exact ground state which is fully phase separated. Thus, we may evaluate the 
ground state energy to find 


Eq — £0 - 

9uv 


N-l 


— VN/2 + ^ 


(S2) 


Requiring that we reproduce the exact ground state energy at (/> = tt thus results in the simple expression for 77 (A): 


V{N) = 


1 

Z^i=l Vi - VN/2 

A3/6-A2/4-A/6’ 


(S3) 


with rji those of Table SI. This procedure yields 77 ( 2 ) = ^l/27r, 77 ( 4 ) = 0.297941, 77 ( 6 ) = 0.246318, and 77 ( 8 ) = 
0.214440. For larger A the exact coefficients are not known, and we instead extrapolate fj towards the thermodynamic 

















2 



FIG. SI. fj{N) as a function of 1/A^. For N = 2,4, 6,8 we use the exact result, indicated by dots, while for larger N we 
approximate it by the dashed line, as outlined in the text [see Eq. (S6)]. The dotted line corresponds to f]{N) calculated from 
the large N limit, Eq. (S5). 


limit N ^ oo. In this limit, the domain wall may be ignored, and the energy to leading non-trivial order is that of a 
Lieb-Liniger gas, Eq. (5) in the main text: 


7V2 128v^ fV5/2 

2 457r2 


(S4) 


This result was obtained in Ref. [42] by applying the local density approximation to the uniform space Bethe ansatz 
result. Hence, in the limit N —>■ oo, the function fj{N) is, to leading order. 


Vll{N) 


64^2 1 


(S5) 


To go to finite iV, we choose the simplest possible extrapolation, requiring that the correction to Eq. (S5) 

produces the correct The result for > 8 is 


fi{N) « 77LL(iV) (1 - 0.063343/iV), 


(S 6 ) 


which is shown in Fig. SI. 


ACCURACY OF THE STRONG-COUPLING ANSATZ 

The high accuracy of the strong-coupling ansatz developed in Ref. [33] and used here to determine the spin-exchange 
coefficients beyond the few-body limit may be confirmed by computing the overlap between the eigenstates |' 0 ) of the 
approximate model and the eigenstates |'0) obtained in the exact solution of the problem. This we have computed for 
balanced systems containing at most 8 particles, analogously to what was done in our earlier work [33]. As may be 
seen in Fig. S2, the wavefunction overlaps for both the ground and first excited state (in the upper branch) remain 
extremely close to 1 for all ratios of intra- and interspecies interaction strengths. The largest deviations (of at most 
a few parts in 10^) are found near the boundaries of the EM region. 


SPIN GAPS 

The stability of the system in the presence of a small inter-conversion coupling may be analyzed in terms of the 
“spin gap”, i.e., the energy difference between the ground state of a system of N particles, and the configuration 
in which one of the spins is flipped. In the EM phase ( 37 r /4 < (f) < 57 r/ 4 ), the ground state corresponds to a fully 
polarized gas of N identical bosons, and we thus define the EM spin gap to be AEfm = Eq{N — 1,1) — i?o(A,0), 
where Eo(A^, is the ground state energy of a gas composed of spin-up particles and A^ spin-down ones. For 





3 



0 

FIG. S2. Wavefunction overlaps between exact and approximate ground (solid) and first excited (dashed) states for = 

2, 3, and 4 particles. 


all other values of 0, including those in the AFM phase, the ground state corresponds to a balanced system with 
= N^. Thus the spin gap in this region is instead defined as AEafm = Eq{N/2 — l,N/2 + 1) — Eo{N/2,N/2). 

At the classical Ising points 0 = 0,7r, we know the ground states exactly and thus we also know the spin gap, 
which at both points is AE = 2fj{N){N — l)/g. This follows from the observation that in both cases the energy 
cost associated with flipping a spin is smallest at the trap edge. We then apply perturbation theory to estimate the 
behaviour of the spin gap around these points. In the AFM regime, we have for |0| <C I, 

AAafm 2 ^ [iV-l-(iV + 3)02] (S7) 

9 

^-^^m(iV)(l-02), (S8) 

where in Eq. (S7) we have kept the leading order correction in N. We display the analytic result, Eq. (S8), in Fig. 3(c) 
to demonstrate the manner in which the AFM spin gap approaches the thermodynamic limit. 

We may perform a similar perturbative expansion around the FM Ising point, but in this case we find that we can 
describe the whole FM region for large N with the following analytic form: 

AEfm ^ 2 (jV- I) ^^^^^\ /cos2,/)-sin^ 0, (S9) 

9 ^ 


as shown in Fig. 3(c). 




















